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Derivation of the procedures that can be apphed in evaluating two-time correlation function 
in terms of coherent-state propagator and corresponding Q-function is presented. On the basis 
that the involved functions are generally exponential in nature, obtaining the two-time second- 
order correlation function is essentially claimed to be reduced to carrying out relatively simple 
integrations. Fundamentally, the time dependence of the operators is transferred to the density 
operator. Moreover, manipulation in reordering the operators is performed by applying the usual 
trace operation. With all details, it is basically observed that the two-time correlation can be readily 
determined once the pertinent coherent-state propagator or Q-function is known. Since working with 
c-number equation is far more handy than the associated operator equation, it is expected that the 
results derived in this contribution can aid in easing the otherwise involving mathematical rigor. 

PACS numbers: 42.50.-p, 42.50.Ar, 42.50.St 



I. INTRODUCTION 

In many conceivable situations, a unique state deter- 
mination through measurement is somewhat physically 
unattainable errand. It must be obvious that detection 
of photon is one of such difficulty. In other words, pho- 
ton measurement may not provide sufficient information 
about the associated state of the photon since the inter- 
action of the radiation with the measuring device is often 
complicated. Even in this context, experimental setups 
for determining the photon number in certain physical 
processes such as atomic excitation is fairly known de- 
spite the compelling fact that the photon number is usu- 
ally very large. It is worth noting that the most com- 
mon method of photon measurement relies on the pho- 
toelectric effect in which the photoelectrons librated by 
the photoionization are counted. It is also good to note 
that the associated detector is sensitive essentially to the 
photoelectrons and fundamentally registers the current or 
the voltage induced by these electrons. As a result, the 
involved devices essentially operate in absorptive mecha- 
nism since the measurement is destructive as the photon 
responsible for the production of the photoelectron dis- 
appears. In addition to this, in the photon detection pro- 
cess, it is imperative assuming that each absorbed photon 
gives rise to no more than one electron and conversely 
each electron is librated only by one photon. 

In actual experimental setup, it is worth noting that 
there can be different ways of measuring photons. One 
of the most obvious procedures in this respect is count- 
ing the number of photons produced by a single source 
with the aid of a single detector. Though this seems 
straightforward approach of measuring the photon, it es- 
sentially lacks a potential for revealing the correlation 
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between the photons produced at different times and po- 
sitions. Remarkably, it is a common knowledge nowa- 
days, similar correlations are responsible for witnessing 
the nonclassical properties of the radiation. In relation 
to this, a number of experiments, basically, interferomet- 
ric by nature have been performed where the photons 
delayed in time and space are counted [l| . Experimental 
arrangements, in which the photons delayed in time are 
counted, are generally designated as delayed coincidence 
measurements. If a light produced by a single source is 
orchestrated to travel over unequal distances like in the 
Michelson-Morley interferometer, an interference pattern 
is produced Sometimes, the photons separated in 
time and space can be counted by more than one detector. 
In such a case, studying the nature of the correlation be- 
tween these photons might be required which makes the 
mechanism of evaluating the correlation between the two 
photons at different time (usually called two-time second- 
order correlation function) an integral part and parcel of 
optical measurements. 

It is, therefore, imperative looking for various ap- 
proaches of mathematically determining the two-time 
second-order correlation function in modern quantum op- 
tics. In light of this, the technique for calculating the two- 
time correlation function based on trajectory approach 
is developed in Q whereas the derivation of its evolu- 
tion equation is presented in [3]. Moreover, the two-time 
correlation function of the harmonic oscillator is evalu- 
ated in 5] . The two-photon quantum correlation between 
Stockes and anti-Stockes radiation in the lambda three- 
level atomic system has also been studied using the mas- 
ter equation formulation and Onsager-Lax regression the- 
orem Furthermore, in recent years theoretical anal- 
ysis of delayed coincidence of the cavity radiation of the 
two-level atom has been addressed 0, H, [1] under vari- 
ous context. Nonclassical features including photon anti- 
bunching and sub-Poisson photon statistics have been re- 
ported upon calculating the two-time second-order cor- 
relation function. There has been also a great deal of 
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interest in studying the statistical and quantum proper- 
ties of the radiation generated by different mechanisms 
in terms of various quantities including second-order [13] 
and multi-time correlations for open and closed systems 

mm- 

In order to extend similar analysis to more involving 
situations, it is presumed to be advantageous introduc- 
ing alternative approaches for evaluating the correlations 
between photons arriving at the site of the detectors at 
separated times. Hence, primarily, the main objective of 
this contribution is directed to fill the gap in this respect. 
To achieve this goal, in the first place, the general con- 
cept of quantum correlation based on the seminal work 
of Glauber [l^] is introduced. In addition to this, the 
procedures of evaluating the two-time correlation func- 
tion applyi ng the most obvious Onsager-Lax theorem 
[TslTlal. the coherent-state propagato r [TtL isl and 

the quasi-statistical distribution functions [11,113, EH, m 

(particularly Q-function) are presented where the former 
is included for the sake of completeness. Taking the al- 
ready available resources in this regard into consideration, 
it is anticipated that the later approaches can be of great 
help in calculating various correlations describablc in the 
form of moments of the radiation at different times. 



II. CORRELATION FUNCTION: GENERAL 
REMARK 

Theory of photon detection requires complete descrip- 
tion of the interaction of radiation with matter, although 
this consideration is not taken into account in this con- 
tribution due to the involved complications and lack of 
complete knowledge of the interaction. In the process 
of photon measurement, since photon absorption mech- 
anism is employed, the detectors are presumed to be in- 
sensitive to the associated spontaneous emission. Conse- 
quently, the annihilation operator of the radiation field 
(E'^'l')) is believed to be the one involved in the counting 
process. In light of this, it is possible to propose that if 
the field undergoes a transition from an initial state (| i)) 
to a final state (| /)) in which a single photon has been 
absorbed, the elements of the transition matrix can take 
the form 
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Assuming the measuring device to be an ideal photon 
detector with frequency independent absorption proba- 
bility, it is not difficult to comprehend that the probabil- 
ity per unit time at which a photon is absorbed at given 
position in space and time can be expressed as [isj 



W,f - (z I ^(-) I /)(/ I ^(t) I 



(2) 



where (E(t))t = 

In actual setting, the final state of the field would not 
be measured. But the measuring device registers the total 



count. In order to obtain the total count, it is necessary 
to sum overall states of the field that can be reached at 
from the initial state via absorption process. In light of 
this, the total counting rate or average field intensity can 
be defined as 

/(r,i)=^(z|^(-)(r,i)|/)(/|^(t)(r,t)|,). (3) 
/ 

With the claim that the final states are complete, / = 
X)/ I/) (/I' the average field intensity can be rewritten as 



/(r,i) = {i I ^(-)(r,i)^(t)(i.,t) I i). 



(4) 



It is not difficult to observe that the expectation value in 
Eq. ([U is quantifiable in view of the initial state alone. 
Moreover, it is noteworthy that in the product of the field 
operators the creation operator precedes the destruction 
operator which corresponds to the normal ordering. 

In principle, it may appear easy and straightforward to 
conceive that recording photon intensities using a single 
detector ensures exhaustive measurement associated with 
the field [H, [l^]. It, hence, turns out to be imperative 
looking for other possible mechanisms for carrying out 
reliable measurement on the field. In this respect, with 
the assumption that there are two fields emerging from 
position r and detected at separate times ti and t2, the 
correlation between the two photons can be quantified 
using the correlation function defined by 

G^^\vMM)=Tr{pE'^-\vM)E'^^\vM)), (5) 

where p is the density operator that describes the state 
of the radiation field and Tr is a shorthand for trace 
operation. Eq. ^ stands for the two-time first-order 
correlation function and it is usually found to be sufficient 
to account for the classical interference experiments. 

It is advisable to resort back to statistical formulation 
since precise knowledge of the field is almost absent. It 
is noticeable that p corresponds to the initial state of the 
radiation field. First and foremost, Eq. ([5]) can be inter- 
preted as the transition probability for the detector atom 
while it absorbs a photon from a field at position r in 
time between t and t + dt. In many instances, stationary 
fields are the common interest in quantum optics whereby 
the correlation function of the field is invariant under the 
displacement of the time variable. Hence the correlation 
function G^^^ (r; ti, ^2) is presumed to depend only on ti 
and t2 through their difference, that is, t = t2 — ti. On 
account of this consideration, it is possible to see that the 
two-time first-order correlation function can be denoted 
as G(i)(r;ti,i2) = G(i)(r;T). 

On the other hand, the joint probability for detecting 
one photoionization at position ri between ti and ti +dti 
and another at r2 between t2 and t2 + dt2 with ti < t2 
is describable by the two-time second-order correlation 
function defined by 
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G(2)(ri,r2;ii,i2) = Tr(p^(-)(ri,ti)^(-)(r2,t2)^(t)(i.2,t2)^(t)(r^^i^)). 

I 



(6) 



It is not difficult to observe that the right hand of Eq. ([6]) 
is time ordered in which the operators at earUer times 
come first and they are also normally ordered wherein 
creation operator comes first. Eq. © generally defines 
the two-time second-order correlation often interpreted 
as the photon delayed coincidences between the two pho- 
tons. 

The discussion up to now is based on the field oper- 
ators. However, in many instances in quantum optics, 
a normalized correlation function in terms of the radia- 
tion or boson operators may be required. In this line, the 
first-order normalized correlation function with the appli- 
cation of the relation between field and boson operators 
turns out to be 



{a'' {t)a{t + t)) 
{af{t)a{t)) • 



(7) 



In the same way, the second-order normalized two-time 
correlation function can be put in terms of creation and 
annihilation boson operators as 



(2). ^ ^ {aHt)a^(,t + T)ait + T)a{t)) 



(8) 



where the radiation field is assumed to be statistically 
stationary. 

It has been a subject of discussion in earlier communi- 
cations that the normalized correlation function is a vital 
tool in identifying the corresponding photon statistics. It 
is common knowledge that when the field satisfies the 
inequality, 



g(2)(r)<.g^^H0) 



(2)^ 



(9) 



for all r less than some critical time r^, the photon ex- 
hibits excess correlation for times less than Tc- This phe- 
nomenon represents photon bunching as the photons tend 
to distribute themselves in bunches rather than at ran- 
dom, since the correlation for photons arriving at the 
same time (r = 0) is greater than the ones coming at 
separated time of r. This actually means that when such 
a light falls on the photon detector more pairs of photons 
are detected closer together than further apart. The re- 
verse of this situation corresponds to the phenomenon of 
photon anti-bunching where fewer photon pairs are de- 
tected closer together. The phenomenon of photon anti- 
bunching, in particular, is one of the possible ways by 
which the nonclassical features of the light is manifested. 
With the aid of this approach, photon anti-bunching is 
found to be a fundamental property that when a two- 
level atom interacts with the radiation in a view that 



a time needs to be elapsed whatever small may be be- 
fore an atom absorbs the radiation after it successfully 
emits it 0, HI]. It is also possible to characterize the 
photon statistics of the light via calculating the two-time 
second-order correlation function. In connection to this, 
it has been known for long that g^^^'ir) = 1 represents 
Poissonian, g'^'^^r) > 1 super-Poissonian and (/'•^■'(r) < 1 
sub-Poissonian photon statistics. 



III. PROCEDURES FOR EVALUATING 
TWO-TIME CORRELATION FUNCTION 

In practice, a solution of the density matrix is not 
always sufficient to determine the two-time correlation 
function. In many instances, it may be required to find 
the transition probability distribution. In some cases, it 
is also possible to evaluate the two-time correlation func- 
tion employing the explicit form of a one-time correlation 
function obtainable with the aid of the master equation 
or the Langevin equation. Although two-time second- 
order correlation function is very important in studying 
statistical and quantum properties of the radiation, its 
evaluation is not often straightforward and easy due to 
the time difference in the operators. Nevertheless, in the 
following, some alternative approaches are discussed in 
order to make the calculation of the two-time second- 
order correlation more easy to handle. 



A. Onsager-Lax regression theorem 

In principle, the correlation function can be readily de- 
rived if the time evolution of the corresponding operator 
is known. This is, essentially, equivalent to the knowl- 
edge of the solution of the Heisenberg or the quantum 
Langevin equation. But, in practice, finding the correla- 
tion of operators evaluated at two different times is not 
obvious as it might appeared from the outset. Quite of- 
ten, in order to calculate the two-time expectation value, 
it is desirable to make use of the utility offered by the 
Onsager-Lax theorem . To this effect, the density oper- 
ator at a time r with r > is expressed in terms of the 
density operator at earlier time i = as 



(10) 



where U{t) is the usual evolution operator defined by 

f7(T) = exp(-i7?s X r), (11) 
in which Hs is the system Hamiltonian. 
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Suppose the system under consideration is consistent 
with Makrovian approximation. This entails that the cor- 
relation between the system and reservoir at equal time 
is unimportant which implies that it is sufficient to write 
Psnit) = Ps{t) ^ PR{t)- In this approximation, the evo- 
lution of a single-time expectation value can be expressed 
following the same reasoning as in Eq. (|10p as 



{A{t + t)) - TrsTrR{UHT)A{t)U(T)ps{t) (g) pR{t))- 

(12) 

Using the cyclic property of trace operation, it is possible 
to see that 

(i(t + r)) = Trs[A{t)TrR{U{T)ps(t) (g) pRmHr))]. 

(13) 

Since Psit)^pRit) represents a single density op- 
erator {psR{t)) that describes the combined system- 
reservoir, it is not difficult to note that 

{Ait + t)) = Trs{A{t)TrR{ps{t + t) (g) pR{t + r))), 

(14) 



which can also be rewritten as 

{A{t + r)) =^G,(r)(i,(i)), 



(15) 



where Gj{Tys are coefficients that depend on r. In the 
same manner, the two-time correlation function can be 
put in the form 

{A{t + T)B{t)) = TrsTrR{U\T)A{t)B{t)U{T) 

^Ps{t)i^PR(t))- (16) 

With the aid of the cyclic property of trace operation, 
one can write 

{A{t + T)B{t)) = Trs(A{t)B{t)TrR{U{T) 

xps{t)(S)PRmHT))). (17) 

Comparing with earlier discussion shows that 

{A{t + T)B(t)) =Y^G,{T){A,(t)B,{t)). (18) 



The procedure of transferring the time r from the oper- 
ator to a function G(t) is described as Onsager-Lax or 
commonly known as quantum regression theorem p^[l6j. 
Fundamentally, what remains is obtaining G(t) based on 
the time evolution of the density operator that depends 
on the underlying physical system and accompanying ex- 
isting circumstance. 



B. coherent-state propagator 

In the previous discussion, the time evolution of the 
quantum system to be studied is described by an opera- 
tor U directly related to the pertinent system Hamilto- 
nian according to Eq. pT|) . On account of the mathe- 
matical difhculty involved in manipulating the operators, 
it is found advantageous employing the corresponding c- 
number equation. One of such formalisms is based on re- 
placing the evolution operator with associated c-number 
function usually designated as coherent-state propagator 
[TtI . [T1 |. With the aid of this approach, theoretical in- 
vestigation of the nonclassical features of the radiation 
generated by parametric oscillator [2^ and spontaneously 
induced entanglement in the cavity radiation of TV two- 
level atoms 27[ have been reported recently. It has been 
observed that this consideration significantly simplifies 
the otherwise cumbersome process. It is, hence, expected 
based on earlier studies that using c-number represen- 
tation can serve the purpose in easing the rigor of deter- 
mining the two-time second-order correlation function. In 
this regard, Shao et al. have used path integral for- 
mulation to find two-time correlation function for systems 
connected with heat bath. With this motivation, in this 
section, the way of obtaining the two-time second-order 
correlation function applying the coherent-state propaga- 
tor would be developed. To this end, an arbitrary func- 
tion (correlation function) of the form 



9{r) 



)a{t)) 



(19) 



is taken for clarity. It is straightforward to realize that 
this expectation value can be expressed in the Heisenberg 
picture in terms of the density operator at initial time 
(p(0)) as 

5(T)=Tr(p(0)at(t + T)a(t)). (20) 

It would not be difficult to notice that the time depen- 
dence can be transferred to the density operator using 
the fact that Tr{p{id)A{t)) = Tr{p{t)A{Q)) as 



g{T)^Tr{p{t)a\T)a) 



(21) 



where p{t) can readily be obtained from /5(0) with the aid 
of Eq. (Iini), that is, 



p{t) = u{t)mu\t). 



(22) 



By introducing the completeness relation for the co- 
herent state, I — J ^\a){a\, in Eq. fH} along with the 
application of Eq. (P^ . it is possible to see that 



9{r) 



(Pa 



Tr{U{t)p{{))U\t)a\T)a\a){a\). (23) 



On the basis of the cyclic property of trace operation, Eq. 
can be rewritten as 



l{r) 



a{a\U{t)p{Q)U\t)a){T)\ 



(24) 
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Assuming the initial state of the system to be described 
by arbitrary state | ao), Eq. ((24|) can be put in the form 

-^a{a\Uit) I ao){ao \ UHt)ix^ {T)\a). (25) 

Now upon defining, K{a,t\(3,0) — {a\U{t)\P), as the 
coherent-state propagator for a single-mode radiation, it 
is found employing the completeness relation for coherent 
state once again that 

-^K{a,t\ai,0){ai\aa). (26) 

Furthermore, with the introduction of the completeness 
relation once again, one can readily see that 

{ao\U\t)a^T)\a) = J {ao\U^ \a2) {a2\aHT)\a) , 

(27) 

which can also be rewritten following the same reasoning 
as 

{ao\UHt)a\T)\a) = f ^^K* ia2,t\m,0) 

J TT TT 

X {ao\a3){a2\a''{T)\a), (28) 



where K*{a2,t\ a3,0) = (ag | UHt) \ aa). 

Following the same line of argument, it is not difficult 
to see that 

{a2\a^{T)\a) ^ Tr{p' {Q)a\T)), (29) 

where p'(0) = |Q)(a2|. Applying the property of trace 
operation, it is possible to shift the time dependence to 
this density operator as we have done before. That is, 

{a2\a\T)\a) ^ Tr{p' [rp), (30) 

where p'(r) — U {T)\a) {a2\U'^ {t) . Inserting the complete- 
ness relation for a coherent state into Eq. ([30]) leads to 

/(P (\ 
—^aX{a4:\p'(T)\aA), (31) 

which can also be rewritten as 

/(P 
aXK{ai, t; a, 0)K*{a4, r; a2, 0). 

(32) 

On account of Eqs. (|25| . (|27| . and (|32| . one readily gets 



9ir) 



(Pa Pai Pa2 (Pa^ Pa^ 



aQ!4(ai|Q!o)(cko|a3)^(a4, T\a, t)K* {a/^,T\a2^Q)K{a, t\ai,Q)K* {a2,t\az,Q)- (33) 



It is not difficult to realize that the correlation function 
can be evaluated using the coherent-state propagator in 
the same manner. What essentially remains to be done is 
to extrapolate this derivation to the case when there are 
four operators instead of two, adapt the coherent-state 
propagator for different variables, and then carry out the 
integration straightaway. It is worth noting that funda- 
mentally similar results have been obtained in [285 ■ Char- 
acteristically, the method of evaluating the coherent-state 
propagator for most general Hamiltonian is provided in 
Ref. [23. It has been observed that the coherent-state 
propagator associated with the quadratic Hamiltonian 
can be expressed in exponential form, which makes the 
involved task quite simple despite the number of integra- 
tions to be performed. For instance, following the pro- 
cedure introduced in [2I], the coherent-state propagator 
for N two-level atom in the cavity and free space [30] and 
parametric oscillation [2^ have heen. calculated and it is 
found to be represented by simple exponential function. 



C. Q- Function 

One of the methods applicable while c-number equa- 
tion is used to study quantum properties instead of the 
pertinent operator equation is the quasi-statistical dis- 
tributions. Quasi-statistical distributions are c-number 
functions related to the density operator in certain pre- 
determined order. One of these functions is the Husimi Q- 
function which corresponds to the normal ordering of the 
density operator. Quite generally, Q-function can be em- 
ployed in calculating various order of moments. In view 
of the earlier efforts, it is expected that making use of the 
advantageous offered by this function can ease the rigor 
of obtaining the two-time second-order correlation func- 
tion. To this effect, it is noticeable that the Q-function 
that corresponds to a time-dependent density operator 
can be defined as 



Q{a) = -{a\U{t)\ao){ao\UHt)\a). 



(34) 
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Introducing the completeness relation for a coherent state 
twice shows that 

If (Pa^ (Puq 
Q[a) = - / 

IT J TT TT 

X {a\U{T)\a:-,){a^\a(3){a(3\aQ){ae\U\t)\a). (35) 

On the basis of the definition of the coherent-state prop- 
agator, one gets straightaway 



1 f (Pa^ (faQ 



TT J TT TT 

xK{a, T\a5,0)K*{a, ijae, 0)(Q!5|ao}(ao|a6}. (36) 

Taking this into account, it is possible to rewrite Eq. ([33)) 
as 



5(r) 



(Pa 



(Pa2(Pa4Q' {a4, a^, T)Q{a, t)aa1. 



(37) 



It is not difficult to realize that the Q-functions in Eq. 
p7p are the pertinent quasi-statical function representing 
the system described in terms of different variables. 

On the other hand, the time evolution of the quantum 
system can be directly obtained from the corresponding 
density operator. In this line, suppose a two-time corre- 
lation function can be expressed as 



g{T) = Tr(a^a{T)p{t))- 



(38) 



It is worth noting that the density operator can be ex- 
panded in the normal order applying the power series 
representation as 



p(<) =^G™(t)at' 



(39) 



Lm 



Therefore, introducing the coherent state completeness 
relation leads to 

5W = / ^5^a„(i)Tr(ata(T) | a){a \ at'a"). 

1,771 

(40) 

In view of the action of the boson operators on the 
coherent state, (a | = a* {a \ and {a \ a'" = 
(a -I- gfr) (a I, it is possible to see that 

X rr(ata(r) I a)(a I). (41) 
On the basis of the fact that 

Q(a,a*,t)-iVa„(t)a*'a", (42) 

TT ^ ^ 



while the operators are initially put in the normal order, 
one finds 

g(r) = J PaQ (^a*,a+^,tj Tr{a^a{T) \ a){a \). 

(43) 

Making use of the cyclic permutation of trace operation 
results in 

Tr{a^a{T) \ a){a \) = a* {a \ a{T) \ a). (44) 

Furthermore, it is possible to express 

{a I a(T) I a) = Tr{a{T)p), (45) 

where p — | a) {a \. With no doubt, the time factor can 
be transferred to the density operator in light of earlier 
discussion. That is, 



{a I a(T) I a) = Tr{ap{T)). 



(46) 



Now upon introducing a coherent state completeness re- 
lation once again, one can write 

{a I a{T) I a) = 1 ^/3Tr(| | p(r)). (47) 

With the aid of the definition of the Q- function in terms 
of the density operator, one can see that 

(a I a(r) | a) = [ d? (3Q{(3, f3\T)(3. (48) 



Hence on account of Eqs. P5|) and one finally ob- 
tains 

(49) 

The Q-functions in Eq. P5|) are also the same Q- 
function pertinent to the quantum system under consid- 
eration in terms of different variables. In the same way 
this approach can be extrapolated when there are more 
than two operators. 



IV. CONCLUSION 

Detailed derivation of various approaches with which 
the two-time correlation function can be evaluated is pre- 
sented. It is assumed that employing c-number equations 
instead of the corresponding operator equations eases the 
involved mathematical rigor. It is basically shown that 
the two-time second-order correlation function can be ob- 
tained from the pertinent coherent-state propagator and 
Q-function. Since the operators, consequently the pho- 
tons, are presumed to be described at different times in 
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present contribution, the coherent-state propagator and 
Q-function representing the quantum system under con- 
sideration should be defined in terms of specially dif- 
ferent time parameters. This entails that the two-time 
second-order correlation function is expressed in terms of 
these functions that can be associated with different al- 
ternatives. In the view that quite significant number of 
quantum systems have quadratic Hamiltonian, the cor- 
responding coherent-state propagator and Q-function are 
claimed to be well behaved exponential functions. There- 
fore, though the number of integrations to be carried out 
are found to be large undoubtedly they reduce to quite 
ordinary standard integrals. The possibility of rewrit- 



ing the approaches following from applying the coherent- 
state propagator in terms of the associated Q-function is 
believed to be essential in verifying the obtained results. 
In the same way, the possibility of rewriting the way of ob- 
taining the two-time second-order correlation making use 
of Q-function in terms of different variables can also pro- 
vide an alternative means of determining it. It is hence 
expected that the detailed derivation presented in this 
work lays a foundation for viable approach of obtaining 
correlations of various moments evaluated at two differ- 
ent times. With no doubt, this procedure can readily be 
employed in evaluating various quantum correlations at 
equal time. 
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